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ABSTRACT 

There  is  a  need  for  a  fast  method  of  analyzing  bathythermograph 
traces  and  this  need  is  approached  by  the  means  of  a  high-speed  digital 
computer.   The  theory  behind  the  computer  program  is  outlined.   Both 
synthetic  and  real  data  cases  are  run  as  examples  using  both  data  card 
decks  and  magnetic  tape  inputs.   The  program  has  been  designed  to  read 
digitized  bathythermograph  traces  and  then  analyze  them  objectively  by 
Gaussian  and  non-Gaussian  methods  for  the  top,  center,  and  base  of  the 
main  thermocline.  Additionally,  such  features  as  multiple  thermoclines, 
inversions,  and  thermal  transients  are  identified  also  and  their  key 
points  are  included  in  the  informational  data  printout. 
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1.    Introduction. 

Since  the  introduction  of  the  bathythermograph  by  Rossby  and 
Montgomery  (1934),  and  its  further  modification  and  perfection  by 
Spilhaus  (1938),  there  has  been  almost  continuous  thermal  structure 
information  funneled  into  the  various  collection  agencies  which  has 
been  handled  manually  under  procedures  outlined  by  LaFond  (1951). 
This  data  collection  has  resulted  in  vast  numbers  of  bathythermograph 
slides  that  have  been  noted,  and  then  stored  away  for  future  analysis. 
In  more  recent  years  bathythermograph  data  has  been  stored  on  aperture 
cards,  using  techniques  similar  to  those  suggested  by  Sauer  (1964). 
The  Canadian  Oceanographic  Data  Center,  Scripps  Institution  of  Oceano- 
graphy, and  the  Fleet  Numerical  Weather  Central,  Monterey,  California, 
have  further  modified  this  procedure  by  a  method  similar  to  that  sug- 
gested by  Sauer  and  Hope  (1967),  where  a  device  has  been  invented  which 
will  take  incoming  bathythermograph  slides  and  convert  them  directly  to 
digitized  data  format.   This  system  described  by  Sauer  and  Hope  (1967) 
consists  of  taking  the  grid/slide  combination  and  transferring  it  to  an 
aperture  card  by  means  of  the  BT  Processor,  a  modified  3M  Processor 
Camera.   It  has  the  advantage  of  preserving  the  original  analog  trace. 
The  aperture  card  is  then  fed  through  a  hopper  device  directly  into  a 
Caps  Jeffree  projector  with  a  D-MAC  curve-follower  unit.   The  time  re- 
quired for  this  unit  to  process  one  bathythermogram  trace  is  approxi- 
mately 45  seconds  with  a  claimed  accuracy  of  1.0%  of  the  full  depth 
range  and  0.1  degrees  Farenheit.  After  exhaustive  testing  at  the 
Naval  Hydrographic  Office  as  described  by  Stewart  (1963)  Sauer  and 
Hope  (1967)  of  the  Canadian  Oceanographic  Data  Center  felt  this  degree 
of  reliability  in  fact  was  not  approachable.   Accordingly,  they 


converted  their  system  to  record  temperature  to  a  minimum  of  every  0.1  de- 
grees Celcius  and  to  the  nearest  whole  meter. 

Papers  by  Tully  (1961,  1964)  show  that  there  are  certain  features 
which  must  be  identified  in  order  to  classify  a  bathythermograph  trace. 
Tully  states  that  these  few  parameters  can  be  used  to  identify  the  thermal 
structure  in  any  area  of  the  ocean.   Among  these  features  is  a  potentially 
isothermal  layer,  which  may  or  may  not  contain  transients.   There  is  a 
seasonal  thermocline,  in  which  the  gradients  may  be  regular  or  irregu- 
lar.  Below  this  layer  is  a  sub -thermocline  layer  in  which  the  gradients 
are  less  than  in  the  thermocline.   Superimposed  on  these  features  are 
various  types  of  short  term  fluctuations  which  possess  varying  degrees 
of  significance,  such  as  temperature  inversions. 

Tully  (1964)  suggests  in  this  same  paper  that  the  thermal  structure 
can  be  defined  by  the  following  features: 

(1)  Surface  Temperature. 

(2)  Potential  Layer  Depth.   This  figure  includes  the  average 
depth  to  the  top  of  the  thermocline,  plus  any  modifications  due  to  inter- 
nal waves,  etc. 

(3)  Depth  to  the  Bottom  of  the  Thermocline. 

In  many  oceanic  areas  this  depth  has  an  associated  distinct  tempera- 
ture minimum  and  therefore  is  very  definitive.   In  areas  like  the 
tropics  it  is  usually  not  a  significant  feature. 

(4)  Thermal  Magnitude  of  Thermocline.  To  researchers  this 
may  be  a  very  significant  feature,  because  it  can  signify  a  boundary 
to  restrict  fish  movement,  water  mass  mixing,  etc. 

(5)  Other  Less  Definable  Features.  There  are  other  parameters 
that  may  be  considered,  like  the  location  of  temperature  minimum  in  a 
thermal  inversion,  mixed  layer  depth,  etc.   These  features  cannot  be 
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fitted  easily  into  a  generalized  machine  processing  format  for  ready 
identification.   They  are  not  common  features,  and  may  have  special 
variables  that  depend  solely  upon  the  area  of  the  ocean  being  investi- 
gated, or  upon  the  season  of  the  year. 

As  Boston  (1966)  states,  "until  the  information  on  bathythermograph 
slides  can  be  conveniently  analyzed  by  machine,  much  information  about 
the  thermal  structure  of  the  upper  layers  of  the  ocean  is  quite  likely 
to  remain  hidden  in  an  ever  growing  mountain  of  smoked  and/or  gold  film 
glass  slides,  photographs  of  same,  or  punched  cards."   Two  distinct 
problems  arise:   (1)  the  digitizing  of  the  bathythermograph  data  to  a 
form  that  can  be  processed  by  high-speed  computers  and  data  handling 
systems;  and  (2)  specification  of  a  program  that  will  take  this  digi- 
tized data  and  process  it  for  use  in  oceanographic  research;  the  pro- 
gram needs  the  generality  to  be  able  to  handle  any  shape  trace,  and  flexi- 
bility so  that  further  routines  may  be  added  or  subtracted  at  will  with- 
out destroying  the  basic  integrity  of  the  program.  As  we  have  seen, 
some  progress  is  being  made  in  the  handling  of  the  first  problem. 

The  second  problem  is  the  topic  of  this  thesis;  the  objective 
analysis  of  digitized  bathy thermogram  traces.   The  standard  features 
as  suggested  by  Tully  (1964)  are  identified  and  printed  out  in  a  concise 
informative  way.  One  or  two  additional  features,  like  the  identifica- 
tion of  thermal  inversions  and  transients,  are  included  in  an  attempt 
to  show  that  further  features  may  be  added  to  this  program  according  to 
the  requirements  of  the  individual  user. 
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2.   Method  of  Analysis  of  Thermal  Structure. 

There  will  be  two  separate  analytical  methods  developed  for  handling 
the  objective  analysis  of  the  thermal  structure.   In  his  technical  paper, 
Boston  (1966)  suggested  the  following  two  methods,  on  which  this  thesis 
is  based.   One  method  consists  of  a  Gaussian  temperature  structure  being 
assumed  for  the  thermal  parameter  discussed.   The  second  method  is  a 
common  finite  difference  method,  utilized  to  find  inflection  points  in 
the  various  trace  slopes. 

A  prerequisite  to  both  methods  is  that  the  data  from  the  original 
analog  bathythermograph  trace  be  reduced  to  an  ordered  set  of  depths 
read  at  some  predetermined  temperature  interval.   From  this  ordered  set, 
the  thermal  structure  can  be  defined  objectively  in  terms  of  certain  pri- 
mary and  secondary  features.   Both  methods  require  only  the  first  and 
second  finite  differences  with  depth;  therefore  only  data  from  depths 
corresponding  to  some  constant  thermal  interval  (conveniently  0.5  to  1.0 
degrees)  are  required  to  handle  the  analytical  technique.   However,  it  is 
convenient  to  include  temperature  data  corresponding  to  these  points,  as 
in  the  computer  output  there  are  locations  where  information  as  to  the 
thermal  magnitudes  are  as  important  as  the  spacial  magnitudes. 

Thermal  structure,  as  defined  in  this  thesis,  is  the  relationship  of 
the  temperature,  T,  to  depth,  Z,  as  a  function  of  time  and  position. 
Furthering  earlier  ideas,  Tully  (1964),  Tully  and  Giovando  (1963),  and 
Uda  (1963)  show  in  their  papers  that  the  thermal  structure  changes  with 
season,  latitude,  and  even  oceanic  region.  However,  Tully  (1964)  also 
suggested  that  the  number  of  different  features  for  a  given  application 
is  small.   Section  A  considers  those  features  common  to  most  oceanic 
regions  within  the  upper  300-400  meters  of  the  ocean.   Section  B  considers 
those  features  that  are  of  major  concern  in  such  fields  as  antisubmarine 
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warfare,  biological  forecasting,  and  fisheries  that  are  not  discussed 
or  included  in  Section  A. 

A.    Primary  Features: 

The  primary  features  as  suggested  by  Tully  (1964)  and  Boston  (1966) 
consist  of  several  convenient  parameters;  (1)  surface  temperature,  (2) 
depth  of  potential  layer  depth  (upper  isothermal),  (3)  depth  of  the 
bottom  of  the  thermocline,  and  (4)  magnitude  of  the  thermocline. 

The  potential  mixed  layer  is  that  as  defined  by  H.  0.  Pub  SP-105, 
Vol  5,  and  is  considered  to  be  isothermal,  although  this  "isothermal" 
region  accepts  small  positive  or  negative  gradients  up  to  0.3  degrees 
Farenheit  per  hundred  feet.   The  termination  of  this  nearly-isothermal 
layer  is  usually  abrupt  enough  to  clearly  mark  the  region  where  the 
thermocline  or  other  feature  takes  place.  The  depth  of  the  mixed  layer 
will  be  defined  as  Zi  the  top  of  the  thermocline.   Likewise,  Zl^l  the 
base  of  the  thermocline,  can  be  considered  the  top  of  the  lower  nearly- 
isothermal  area. 

Once  the  top  and  bottom  of  the  thermocline  have  been  defined,  with 
It   as  the  temperature  at  2-j.  ,  and  I  b£_as  the  temperature  at  ^^V  » 
the  definition  of  the  magnitude  of  the  thermocline  is 


AT  =  Tt-Tbt 


(i) 


B.    Secondary  Features 

The  secondary  features  are  usually  considered  as  time  varitions  in 

the  primary  features  of  Section  A,  and  can  occur  at  any  point  in  the 

vertical  water  column  being  considered. 
1.    Transients 

The  term  transient  will  be  defined  here  as  referring  to  short  term 
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fluctuations  of  the  thermal  structure.   It  is  usually  the  result  of  wind 
mixing,  and  vanishes  over  a  period  of  a  day  or  more.   The  program  pro- 
duced for  this  thesis  will  recognize  these  transients,  assuming  that 
these  thermal  transients  take  the  usual  form  of  small  but  sharp  gradients 
that  appear  as  irregularities  on  the  main  thermal  structure.   Transients 
are  often  found  in  association  with  the  well  known  "afternoon  effect." 
This  program  will  aid  in  the  verification  of  the  prediction  methods  out- 
lined by  Tabata  and  Giovando  (1962). 

2.  Multiple  Thermoclines 

Heating  and  wind  mixing  conditions  can  also  produce  more  than  one 
stable  thermocline.   The  growth  of  two  stable  thermoclines  has  been  well 
documented  by  Tully  and  Giovando  (1963).  The  occurrence  of  these  multiple 
thermoclines  is  sufficiently  common  and  vital  to  underwater  acoustic  pre- 
diction systems  that  they  merit  identification. 

3.  Inversions 

Inversions  are  a  common  oceanic  feature,  an  example  of  which  has 
been  described  by  Roden  (1964),  and  must  be  identified  if  one  is  to  ac- 
curately predict  underwater  sound  transmission.   The  thermal  inversion 
takes  form  when  the  temperature  below  the  main  thermocline  (below  Zfc>t) 
decreases  with  depth  to  a  minimum,  then  increases  markedly,  and  finally 
decreases  again,  but  now  usually  at  a  more  reduced  rate  than  the  first 
decrease.  These  thermal  features  are  of  extreme  interest  in  underwater 
communications  as  the  point  of  minimum  temperature  in  this  thermal 
structure  forms  the  axis  of  a  "sound  channel"  formed  by  the  thermal  bend- 
ing of  acoustic  rays  by  this  inversion. 

The  following  methods  will  be  used  to  identify  these  primary  and 
secondary  features  within  the  computer  program  designed  for  this  thesis. 
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C.    The  Gaussian  Thermocline 

A  gaussian,  or  normal,  distribution  (see  Figure  1)  of  temperature, 
T,  as  a  function  of  depth,  Z,  is  given  by 

T(Z)  =  _L_    e     dZ  (2) 


whose  corresponding  frequency  function  is 

dial  z  4=.  €  <3> 

If  the  distribution  function  is  T \cL~ V(\\    ,  where  G~     greater  than  0  and 

V  cr    J 
m  are  constants,  then  the  frequency  function  becomes 

djT(z)    -    \    -  e    -**■  (4) 

dZ     crfi™ 

The  function  described  by  (4)  has  the  following  characteristics: 

(a)  it  is  symmetric  about  Z  ■  m, 

(b)  it  has  two  symmetric  points  of  inflection  at 
'Z—ta*.  &    ,  and 

(c)  it  has  a  maximum  rate  of  change  at  2  s^!Cr|5 

A  change  in  the  numerical  value  of  m  causes  a  displacement  of  the  curve 
in  the  vertical  direction,  Z,  but  does  not  alter  the  form  of  the  curve. 
A  change  in  O"  has  the  effect  of  altering  the  scale  of  both  coordinate 
axes.   The  smaller  the  value  of  CT     ,  the  more  concentrated  the  curve  is 
about  the  point  Z  =  m„   If  the  water  surface  is  defined  as  Z  =  0,  and  the 
depth  considered  to  be  infinite,  then 


r°     „__a 


■rta)-  J. 


ffHf  / 


-(r-»r  (5) 


■00 

is  the  function  defining  the  vertical  temperature  distribution.   The 
center  of  the  thermocline,  or  the  center  of  the  curve  of  the  llfcJjQOj 

distribution,  is  at 
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yr '       t  i        "t  1  1  l  i 

Normalized  Temperature 


dZ 


-m+J3< 


-m 


-m-JScr 


L¥  I 


i       i       i     0*5     i       i       i       i LQj 

(after  Boston  1966) 

Fig.1     Gaussian  thermocline  and  its  frequency 

function  (first  derivative) 


i<> 


Z  =ZC  =  -w  (6) 

The  top  and  bottom  of  the  thermocline  are  then  defined  as  being  located 
at  the  points  of  the  maximum  rate  of  change  of  slope  of  the  T(  aLxJ22  J  curve; 
namely, 

Zt  =  -m  +  a-JT  (7) 

^bt*  "^  ~  ^^  (8) 

which  are  illustrated  in  Figure  1. 

The  handling  of  bathythermograph  data  leading  to  the  application  of 
results  (6),  (7),  and  (8)  is  very  simple,  since  values  of  m  and  O"  may  be 
computed  directly  from  the  (t,Z)  distribution  curve.   Consider  the  first 
two  moments  of  the  (T,Z)  curve  about  the  origin,  Z=0,  whereby  the  (nth) 
moment  of  the  curve  about  the  origin  is  meant  the  mean  values  of  Z^ 
defined  as  Z     >  where  n  is  a  positive  integer.   Then  with  the  normal 
curve,  the  first  moment  is  identical  with  the  center  of  the  thermocline: 


2c  = 


ZCrt  dT  (9) 


Ttol  -T(-oo)  kk») 

where  T(0)  is  the  temperature  at  the  surface  Z  =  0  and  T  (-oo)  is  the 
temperature  at  the  bottom  of  the  lower  isothermal  layer  which,  in  this 
idealized  case,  is  assumed  to  extend  to  the  bottom  of  the  ocean.   The 
second  moment 


z*« 


Tkrt-Tf-coUrt.^ 


ZzclT  do) 


2 
can  be  used  to  compute  CT     ,  since  the  variance,  0"  »  is  given  by 
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In  practice,  several  simplifications  are  possible.   The  temperature  T  (-co) 
is  measured  at  the  bottom  of  the  bathythermograph  cast.   The  differentials, 
dT,  become  finite  differences,  &Tt   and  the  integrals  are  replaced  by  sum- 
mations.  If  the  depths  of  the  temperatures  are  read  at  one  degree  intervals, 
for  example,  then  the  T's  can  be  replaced  by  unity. 

One  must  remember  that  although  integration  is  a  continuous  process, 
the  less  accurate  summation  over  finite  intervals  is  not.   The  larger  the 
interval,  the  less  accurate  is  the  summation.   The  working  value  of  sigma, 
as  computed  by  the  routine,  cannot  exactly  equal  the  correct  value  obtained 
from  integration  over  the  normal  curve.   Thus  equations  (7)  and  (8)  must 
be  modified  by  some  weighing  factor  depending  upon  the  temperature  inter- 
val chosen.  Allowing  this  factor  to  be  defined  as,  sk,  equations  (7)  and 
(8)  become 

2.^.    -  -nr\  +  skjar  (7A) 

2bt  ~  -m  -  skJa  (8A) 

where  a  temperature  interval  of  0.5  degrees  gives   sk  equal  to  1.47  and  a 
temperature  interval  equal  to  1.0  degrees  gives  a  value  for  sk  equal  to 
1.30. 

In  summary,  the  objective  definition  of  the  top,  bottom  and  center 
of  a  Gaussian  thermocline  requires  only  readings  of  (a)  surface  temperature, 
(b)  bottom  temperature,  and  (c)  depths  at  some  constant  pre-chosen  tempera- 
ture interval. 

The  center  of  the  Gaussian  thermocline  is 

\?-,,T\        (12) 


zt-  zcrjj.  *  —4- — )  z(t„v 

1 1  -  I  bt  (— » 
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Since 

-,  i-P-  \  \  N"1 v    2 

(13) 


z^-T^-.EzrrnW 


i  -    'bt 

Then 


cr*  -  2TU1  -  Zc1 


(14) 


The  top  and  bottom  of  the  Gaussian  thermocline  readily  follow  from  (7A) 
and  (8A). 

One  problem  remains"   the  decision  as  to  whether  or  not  a  normal 
distribution  exists;  lacking  this  distribution,  the  inaccuracies  built 
into  this  process  are  not  tolerable.   This  determination  can  be  accomp- 
lished by  computation  of  the  third  moment.   However,  in  the  following 
section  it  is  apparent  that  this  computation  is  not  necessary. 

D.    The  non-gaussian  Thermocline 

The  analysis  of  the  more  common  type  of  thermal  structure,  the  non- 
gaussian  thermocline,  requires  the  computation  of  finite  differences  on 
Z(Tn).   The  first  finite  difference  will  be  defined  as 

taCTn)  -  ZCO  -Z(Tn)  (15) 

and  the  second  finite  difference  as 

A* 2(Tn)  '  b&Tw)  -  oZiTA  (i6) 

Generally,  the  first  finite  difference  is  always  positive.   The  second 
finite  difference  will  be  negative  until,  in  this  case,  the  point  of 
inflection  of  the  (T,Z)  distribution  curve  is  reached,  at  which  point 

&2I(Tf\j      becomes  non-negative.      If   the   first  non-negative   second 
finite   difference    is    AZ(Tio)       s    then  the   center  of   the    thermocline   is 
ZlTn^        . 
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In  order  to  locate  Z  ,  only  the  upper  portion  of  the  (T,Z)  distribu- 
tion curve  is  considered,  where  a  (£_  —  Zc.  J  synthetic  lower  portion  to 
this  curve  segment  is  constructed  so  that  together  with  the  upper  real 
segment  a  normal  distribution  is  formed.   From  this  curve,  0"   is  com- 
puted by  methods  described  in  the  previous  section  and  Z   is  given  by 
formula  (7A) .   (see  Figure  2b). 

Similarly,  only  the  lower  (Z— "Zc^  part  of  the  (T,Z)  distribution 
curve  is  utilized  in  computing  Z ^  .   An  upper  portion  of  the  curve 
is  added  to  the  segment  and  a  new  value  of  C"      is  computed.   ^bL^s 
then  calculated  by  using  formula  (8A) .   (See  Figure  2c). 

E.    Irregular  Thermal  Features. 
1.    Transients 

Transients,  according  to  Boston  (1966),  are  usually  considered 
to  be  less  than  one  degree  in  magnitude  and  might  not  be  seen  by  a  greater 
thermal  interval.   Thus,  it  is  necessary  to  have  a  small  grid  spacing, 
say  0.1  degrees,  to  identify  these  transients.   Spilhaus  (1938)  relates 
that  the  bathythermograph  is  seldom  more  accurate  than  0.5  degrees, 
borne  out  by  the  testing  of  Stewart  (1963),  due  to  such  problems  as 
dull  stylus,  hysteresis,  and  calibration.   Even  if  the  instrument  were 
exact,  and  hysteresis  were  not  present,  the  thickness  of  the  trace  limits 
reading  accuracy  to  about  0.2  degrees.   Hence,  although  the  grid  spacing 
chosen  for  this  program  is  unrealistic  at  0.1  degrees,  it  is  retained  with 
the  possibility  of  more  accurate  thermal  instruments,  or  more  accurate 
bathythermograph  readings,  in  the  future. 

If  such  accuracy  were  possible  (see  Figure  3)  transients  could  be 
identified  by  simply  identifying  the  sign  changes  of  A"Z.(Trt)  -   If  the 
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Fig.  3      Normal  and  expanded  scale  transient 
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first  negative  value  encountered  was  &Z.('n)  then  the  beginning  of  the 
first  transient  would  be  at   Z(T/\*i)  .   Then  the  end  of  the  first  thermal 
transient  and  the  beginning  of  the  second  transient  is  signified  by  a 
sign  change  of  AZ(TVi)  from  positive  to  negative.   If  this  occurs  at 

A~ZOV>;  then  the  end  of  the  first  transient  occurs  at   Z(Tn^    and 
the  beginning  of  the  next  occurs  at  Z.(T/Ht)   •   This  process  will  con- 
tinue until  the  top  of  the  first  major  thermocline  is  encountered; 
generally  we  are  interested  in  only  the  transients  in  the  upper  isothermal 
region. 

2.   Multiple  thermoclines 

The  analysis  of  a  thermal  distribution  containing  more  than  one 
major  thermocline  proceeds  along  the  lines  just  discussed  in  the  pre- 
vious section  on  transients.   In  Boston's  article  (1966)  he  states  that 
a  grid  spacing  of  one  degree  is  generally  adequate  if  there  are  only  two 
thermoclines,  and  a  grid  spacing  of  0.5  is  necessary  should  there  be  more 
thermoclines  present. 

As  in  the  transient  analysis  the  sign  changes  of  AXCTVi)  are 
examined.   If  the  first  positive  value  is   &"2(TV»)   then  the  center  of 
the  first  thermocline  would  be  located  at  Z  (T^,^  )  =Zci  •   Tne  center  of 
the  second  thermocline,  Zc2.   »  ^s  indicated  when  the  sign  of  &"ZjCTV\^ 
changes  from  negative  to  positive.   If  this  occurs  at  £i"2.C"TV>)  >  then, 
as  in  the  first  thermocline,  the  center  is  located  at  "Z (Tn* \)  ~  "2.Ci. • 
The  centers  of  all  subsequent  thermoclines  are  found  by  continuing  this 
process  throughout  the  length  of  the  trace. 

When  the  centers  of  the  shallowest  and  deepest  thermoclines  have  been 
identified,  the  method  outlined  in  Section  D,  the  non-gaussian  thermocline 
will  be  used  to  determine  the  top  of  the  shallowest  thermocline  and  the 
bottom  of  the  deepest  thermocline.   Then  considering  the  fact  that  &"ZOV\) 
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is  a  measure  of  the  slope  of  the  (TSZ)  curve,  the  maxima  in  AZ(Tn) 
corresponds  to  the  most  nearly  isothermal  conditions.   Thus  these  maxima 
can  be  used  to  separate  the  thermoclines.  Once  this  seperation  has  been  ac- 
complished, the  procedure  from  Section  D  can  be  used  again  to  locate  the 
top,  bottom,  and  center  of  each  individual  thermocline.   Considering  the 
fact  that  the  thermal  structure  is  assumed  to  be  continuous,  with  only 
small  spacing  between  consecutive  thermoclines,  there  is  little  need  in 
making  the  identification  any  more  accurate,  or  finer,  than  the  subdivi- 
sions that  are  indicated  by  the  maxima.   Figure  4  shows  an  example  of 
real  traces  in  the  ocean  with  this  multiple  thermocline  feature. 

3.    Inversions 

Large  inversions  (+  0.5  degrees  Celcius)  are  usually  found 
where  water  sources  of  different  origin  are  close  together  or  where 
strong  velocity  shears  are  encountered.   Tully  and  Giovando  (1963)  docu- 
ment this  feature.   Uda  (1963)  and  Tully  and  Dodimead  (1963)  also  show 
that  in  regions  such  as  the  subarctic,  where  there  is  excess  precipita- 
tion, the  summer  profile  may  be  characterized  by  a  nearly  permanent  sum- 
mer inversion  located  at  about  100  meters  depth. 

The  main  feature  (see  Figure  5)  of  an  inversion  is  the  sign 
change  of  the  slope  of  the  (T,Z)  distribution  curve.   The  major  points  of 
interest  to  the  investigator  are  usually  considered  to  be: 

(a)  upper  and  lower  boundries  denoted  by  Zci  and  Zc2 
respectively, 

(b)  magnitude  of  the  inversion,  ZUPP  and  TUPP,  being  the 
spacial  and  thermal  magnitudes,  respectively, 

(c)  minimum  temperature,  TMIN, 

(d)  depth  of  minimum  temperature  (depth  of  axis  of  sound 

channel),  ZMIN. 
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Fig.  5        Typical  BT  trace  with  an  Inversion 
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Referring  to  Figure  5,  it  should  be  noted  that  this  system  of  analysis  will 

require  a  more  elaborate  notation  system,  since  the  curve  is  no  longer 

single  valued  (that  is,  more  than  one  depth,  Z  ,  will  have  a  specified 

n 

temperature) .   In  order  to  distinguish  between  these  depths  the  notation 

Zj  (T~r\\      will  be  introduced  where  the  subscript  j  refers  to  the  order 
in  which  a  multiple  occurrence  happens  from  the  surface  down.   Thus,  if  Tn 
occurs  at  three  depths,  the  shallowest  depth  would  be  indicated  asZ((To)» 
the  intermediate  depth  Z^Cn^  »  an^  tne  deepest  as   Zt,(TVi)  .   As  the 
trace  is  read  vertically  from  the  surface  down,  the  trace  appears  to  be 
single-valued  (one  depth  has  one  temperature),  then  double-valued  (one 
temperature  at  three  depths),  before  returning  to  a  single-valued  state 
again.   Since  the  finite  differences  must  be  computed,  its  first  and 
second  finite  differences  will  be  defined  as  kZj(T"n]  and  A  ZjCTV>), 
respectively. 

The  upper  and  lower  boundaries  of  the  inversion  occur  at  the  first 
double-valued  temperature.   If  this  temperature  is  Tn,  then  Zci  =  ZjOV\) 
an(*  Zc<l  =  Z  2.^01  ■   The  last  double-valued  temperature,  or  a  lower- 
valued  single  temperature  included  between  the  last  double-valued  tempera- 
tures, is  the  minimum  temperature,  TMIN,  and  the  depth  of  this  minimum 
temperature,  ZMIN,  is  defined  from  Z.lTwurt^  •   The  magnitude  of  the 
thermal  inversion,  TUPP,  is  given  by  the  absolute  value  of  the  difference 
between  the  first  and  last  double-valued  temperatures.   Once  these  in- 
version features  have  been  identified,  pertinent  features  above  and  be- 
low the  inversion  still  must  be  identified  by  the  methods  outlined  in 
Section  D. 
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3.    Example  of  Input 

Two  methods  of  input  to  this  program  are  data  deck  and  tape  deck. 
These  two  examples  illustrate  these  methods. 
A.   Data  Deck  Input 

In  using  a  data  deck  input,  where  an  analyzer  does  not  have 
access  to  the  digitizing  machines  and  must  key-punch  his  data  deck,  it 
is  possible  to  input  both  positional  and  curve  data.   The  data  deck  must 
be  made  up  as  follows  (see  Figure  6); 

(1)  The  data  must  be  printed  from  the  bathythermograph 
curve,  at  some  pre-chosen  thermal  interval,  i.e.,  -  0.5  degrees.   The 
data  are  punched  on  the  data  aperture  cards  by  punching  the  temperature, 
followed  by  its  depth  as  picked  off  at  the  chosen  delta-T.   Each  data 
bit  is  allowed  ten  spaces  on  the  aperture  card  and  no  other  data  can  be 
included  in  this  ten  space  unit,  so  that  a  maximum  of  four  data  sets 

can  be  punched  per  data  card,  (since  there  are  only  80  spaces  allowed  per 
card) .   The  thermal  data  is  continued  through  as  many  aperture  cards  as 
necessary  to  insure  complete  recording  of  the  entire  trace. 

(2)  Immediately  following  the  completion  of  the  first 
trace,  the  operator  must  insert  "Control  Card  1",  which  is  a  dummy  con- 
trol card  consisting  of  "0.0"  substituted  for  the  temperature,  T,  and 
any  "non-zero"  value  substituted  for  the  depth,  Z.   This  dummy  card 
signals  the  main  program  that  the  reading  routine  has  completed  reading 
the  data  trace  and  is  ready  to  read  the  positional  data,  if  it  exists, 
or  proceed  on  with  the  analysis  of  the  trace. 

(3)  Following  "Control  Card  1"  we  insert  the  positional 
data  card,  if  we  desire  to  read  in  this  information.   This  card  contains 
the  geographical  position,  both  latitude  and  longitude  to  the  nearest 
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hundredth  of  a  minute,  the  local  date,  the  Greenwhich  Mean  Time  and  date 
of  the  bathythermograph  cast,  and  the  serial  number  of  the  bathythermo- 
graph.  If  the  positional  data  card  is  used,  then  every  position  must  be 
filled  in  exactly  as  noted  in  Figure  6  and  Appendix  II.   If  the  neces- 
sary data  is  not  known,  then  the  unknown  spaces  must  be  filled  with  M0'"s. 
The  format  on  this  card  is  absolutely  vital  and  inflexible. 

(4)  Next  we  repeat  the  above  procedures  (1)  through  (3) 
for  each  bathythermograph  trace  desired.   It  is  not  necessary  to  insert 
a  positional  card  with  each  trace,  but  if  the  positional  card  is  used  it 
must  be  complete,  as  mentioned  above. 

(5)  When  the  entire  data  package  has  been  made  up  of  trace 
data  and  positional  data  cards,  then  "Control  Card  2"  is  inserted  as  the 
last  card  to  signal  the  program  that  the  data  is  finished  and  the  pro- 
gram can  terminate.   "Control  Card  2"  consists  of  a  normal  data  card  for- 
mat, except  the  "0.0"  is  inserted  in  position  for  both  temperature,  T, 
and  depth,  Z.   This  is  always  the  last  card  in  the  data  deck  (if  we  do 
not  consider  the  programming  control  cards  used  by  the  system),  and  will 
always  follow  the  positional  data  card  for  the  last  trace  or  if  this  data 
is  not  used,  then  it  will  follow  the  last  "Control  Card  1."  The  entire 
data  deck  is  assembled  as  shown  in  Figure  6,  and  may  be  of  any  length. 
The  only  programming  prerequisites  are  that  the  thermal  inverval  be  a 
constant  for  each  individual  trace,  although  it  is  not  necessary  for  all 
the  traces  to  have  the  same  interval,  and  the  input  format,  Figure  6, 
must  be  followed  rigorously.   All  other  necessary  computations,  such  as 
sigma,  weighing  factor,  etc.,  are  contained  or  computed  within  the 
program. 


30 


B.    Magnetic  Tape  Data 

When  the  input  data  can  be  digitized,  by  such  facilities  as  the  Fleet 
Numerical  Weather  Central,  Monterey,  or  the  Canadian  Oceanographic  Data 
Center,  a  few  modifications  must  be  made  to  the  basic  program  reading  meth- 
ods as  they  presently  exist  in  Appendix  II.  In  order  to  use  magnetic  tape 
input  it  is  necessary  to  remove  the  cards  that  are  serialized  (right  hand 
margin)  0185  -  0204  and  insert  in  their  place  new  cards  designed  to  read 
the  magnetic  tape.   It  is  important  to  note  here  that  the  format  of  the 
tapes  and  even  the  type  of  tapes  makes  a  vast  difference  in  the  make  up 
of  the  reading  cards  and  great  care  must  be  used  in  inserting  the  cor- 
rect control  and  reading  cards  in  order  to  insure  not  shutting  down  the 
computer.   Before  using  tapes,  be  sure  to  check  the  individual  instruc- 
tion manual  and  the  individual  tape  specifications  and  rigorously  follow 
the  instructions  on  the  format  for  the  make  up  of  these  new  cards  to  be 
inserted.   The  following  information  must  be  considered  in  each  case: 

(a)  whether  seven  or  nine  track  tape, 

(b)  data  density, 

(c)  whether  even  or  odd  parity, 

(d)  whether  the  data  are  translated  or  not, 

(e)  if  translated,  what  machine  language, 

(f)  record  length, 

(g)  whether  or  not  the  data  are  blocked, 
(h)  if  blocked,  what  size. 

All  this  information  is  vital  in  making  up  the  correct  system  control 
cards  for  the  tape  reading  unit.  Once  the  system  control  cards  are  cor- 
rectly made  up  and  inserted  around  the  program,  then  a  normal  "read"  card, 
for  example; 


31 


READ(8,1,END=00) 
1    FORMAT (8A5) 

is  used,  where,  in  this  example,  we  are  reading  on  tape  reader  "8", 

sixty  traces  consisting  of  eight  data  sets  per  trace. 
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4.    Computer  Program  Operation, 

The  computer  first  reads  in  the  data  for  an  individual  trace,  as 
outlined  in  Section  3,  and  then  reads  in  the  positional  data  (if  existent 
in  data  deck) .   Once  the  individual  trace  is  read  into  the  program,  the 
computer  takes  over,  using  the  analysis  method  programmed  by  this  routine 
to  start  the  investigations  for  the  features  outlined  in  Section  2.   Be- 
fore starting  any  further  investigations  the  program  looks  at  the  thermal 
interval  and  choses  one  of  three  possible  values  for  the  weighing  factor, 
SK.   This  value  is  then  carried  throughout  the  calculations  for  this  trace. 
If  the  positional  data  is  included,  the  first  output  is  a  formalized  print- 
out of  this  informational  heading,,   Once  this  has  been  executed,  the  main 
program  starts  a  searching  routine  to  identify  various  possible  features 
before  calling  the  specific  sub-routines  necessary  to  analyze  this  data. 
The  program  first  looks  for  double-valued  temperatures  in  the  vertical 
dimension.   If  these  double-values  exist  the  program  calls  the  inversion 
subroutine,  INVERS,  to  calculate  the  top  and  bottom  of  the  inversion  and 
its  spacial  and  thermal  magnitude.   This  subroutine  first  calls  a  special- 
ized subroutine,  0RDER1,  that  takes  all  the  input  data  from  that  trace 
and  checks  to  insure  that  the  data  is  in  correct  order  from  the  surface 
down  (consecutive  readings)  and  if  it  is  not,  reorders  the  data  in  the 
correct  sequence.   Once  the  data  is  in  correct  order,  then  the  program 
does  a  quick  safety  check  to  insure  that  it  does  not  have  a  non- inversion, 
(where  the  temperatures  first  increases  with  depth  and  then  decreases, 
giving  not  an  inversion,  but  a  positive  surface  gradient  followed  by  the 
normal  main  thermocline,  a  feature  that  the  main  program  might  identify 
as  an  inversion).   If  this  condition  does  exist,  the  subroutine  terminates 
and  returns  the  flow  to  the  main  program  for  further  investigations.   If 
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an  inversion  exists  ,then,  by  the  use  of  do-loops ,  INVERS  calculates  the 
top,  bottom,  and  magnitudes  of  the  inversion  as  outlined  in  Section  E.3. 
Once  the  data  has  been  analyzed  and  these  results  are  available,  the 
subroutine  then  calls  four  specialized  printing  subroutines,  PRINT3 
through  PRINT6,  to  print  out  the  data  in  a  useful  format.   INVERS  then 
returns  the  order  flow  to  the  MAIN  program  for  further  investigations. 

The  MAIN  program  then  investigates  the  thermal  interval,  delta-T, 
to  insure  that  if  the  interval  is  less  than  0.2  degrees  a  thorough 
investigation  for  the  presence  of  transients  is  accomplished.   If  the 
interval  is  0.2  or  less,  the  MAIN  routine  calls  the  transient  investiga- 
tion routine,  TRANS,  which  proceeds  to  analyze  for  the  presence  of  transi- 
ents as  outlined  in  E.l.   Upon  completing  the  analysis,  the  results  are 
printed  out,  if  any  exist,  and  the  flow  returns  to  the  MAIN  routine. 

Once  returned  to  the  MAIN  program,  the  final  investigation  is  made 
into  the  possibility  of  multiple  the rmoc lines.   The  MAIN  program  calls 
the  multiple  thermocline  subroutine,  MANYT,  which  utilizes  analysis 
procedures  outlined  in  Section  E.2.   Within  this  subroutine  is  the  pro- 
vision to  go  to  another  subroutine,  NONGA,  designed  to  handle  the  more 
specific  case,  should  only  a  single  main  thermocline  exist.   The  NONGA 
subroutine  uses  the  procedure  outlined  in  Section  D  to  find  the  top, 
center,  and  bottom  of  the  thermocline.   In  both  subroutines,  looking  at 
thermoclines,  the  results  are  printed  out  by  utilizing  the  specialized 
printing  subroutines,  PRINT1  and  PRINT2. 

A  flow  chart  (simplified)  for  the  entire  process  is  shown  in  Ap- 
pendix I,  and  the  entire  program  can  be  inspected  by  referring  to  Ap- 
pendix II.   It  is  worthy  of  note,  that  this  program  will  not  identify  all 
three  features,  transients,  multiple  (or  single)  thermoclines,  and  in- 
versions, in  a  single  trace  as  the  existance  of  all  three  would  be  very 
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rare.   However,  the  program  will  look  for  transients  in  the  upper  near- 
isothermal  layer  and  then  investigate  the  possibility  of  an  inversion  be- 
low this  area,  or  if  the  inversion  does  not  exist  then  look  at  the  thermo- 
cline(s)  in  the  region  below  the  near  isothermal  layer. 

Also,  as  an  additional  feature  the  program  can  be  called  to  print 
out  a  graphic  display  of  the  input  data  (Appendix  III),  thus  giving  a 
visual  representation  of  the  data  that  is  printed  out  on  the  output  data 
sheets.   This  auxiliary  feature  can  be  called  by  inserting  the  following 
cards  in  the  MAIN  program  immediately  following  the  card  serialized 0199. 

CALL  SEE13(T9Z9ICOUNT) 
which  will  call  the  SEE13  subroutine  that  utilizes  the  library  pro- 
gram, "DRAW",  to  give  graphic  display  of  the  input  data  set  and  fits  a 
curve  to  these  points.   Before  using  this  call,  the  user  must  exercise 
care  and  carefully  read  the  rules  for  usage  of  "DRAW"  in  subroutine 
"SEE13"  as  are  outlined  in  Appendix  III. 

Another  feature  included  as  part  of  this  thesis  package,  is  a  pro- 
visional set  of  cards  (Appendix  IV)  that  will  convert  the  program  as 
contained  in  Appendix  II  into  the  slightly  different  machine  language 
required  when  the  program  is  used  on  the  CDC  6500  model  computers,  in- 
stead of  the  IBM  630  models  for  which  this  program  is  written.   Before 
using  this  conversion  program  two  steps  must  be  taken, 

(1)  all  data  and  system  control  cards  must  be  removed 
and  the  remaining  portions  of  the  included  program 
become  the  data  deck  for  "BCD" 

(2)  All  output  formats  that  involve  printed  words  or 
combinations  of  printed  words  and  numbers  must  be 
converted  to  "H"  format. 
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5.    Output  from  the  Computer. 

The  output  from  the  computer  consists  of  several  possible  combina- 
tions of  data.   Included  are  examples  of  a  few  of  the  possible  combina- 
tions.  All  the  output  from  a  given  trace  is  printed  on  a  single  indivi- 
dual page,  headed  by  the  trace  number  within  the  data  sequence  and  the 
positional  data  for  that  slide  (if  any).   Figure  7  shows  one  of  the  sets 
of  real  data  that  was  taken  from  some  data  traces  made  by  Tully  in  the 
vicinity  of  Ocean  Station  "Papa"  in  1963.   Figure  7  is  a  hand-drawn  en- 
largement of  a  real  test  trace  containing  both  transients  and  a  single 
main  thermocline,  typical  of  that  region,  used  in  testing  out  this  pro- 
gram.  Figure  8  is  a  typical  output  from  the  program's  analysis  of  this 
trace  and  lists  the  positional  data  (which  was  partially  available)  and 
the  date  and  time  of  the  trace.   The  program  then  reveals  the  existence 
of  two  transients,  by  printing  out  the  location  of  the  start  and  end  of 
each  transient.   It  also  gives  the  information  on  the  single  major  thermo- 
cline.  In  this  case  the  data  was  read  in  from  a  data  dekc  after  being 
transcribed  from  the  trace  onto  aperture  cards.   The  data  interval  was 
taken  in  this  case  to  be  0.1  degrees.   Figure  9  shows  the  typical  print- 
out of  a  hand-constructed  case,  also  used  to  test  program,  centered 
around  a  thermal  inversion  (Boston  1966).   Figure  10  is  the  data  output 
for  this  created  case.   Figure  12  is  an  example  of  the  printed  output 
that  would  be  expected  from  a  trace  (see  Figure  11)  containing  multiple 
thermoclines.   In  this  case  two  major  thermoclines  exist  in  the  trace. 
For  interest,  Figure  13  shows  the  output  from  the  drawing  routine  called 
from  SEE13.   It  can  be  seen  that  computer  drawn  trace  of  Figure  13  is  an 
adequate  enough  representation  of  Figure  12  to  enable  the  user  to  use 
this  type  of  computer  output  as  a  rough  aid  in  insuring  correct  working 
of  the  computer  program  or  other  rough  calculations  where  a  visual  presenta- 
tion might  be  desirable. 
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Fig.  7    Ocean  Station  'Papa'     — 

23  April  1963 


(after  Tully  1963) 
BT  Trace 
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START  OF  INVESTIGATIONS  IN  BATHYTHERMOGRAPH  TRACE  NUMBER  4 

GEOGRAPHIC  POSITION  OF  BATHYTHERMOGRAPH  TRACE: 

LATITUDE   .........  36-36.1  NORTH 

LONGITUDE    113-44.7  WEST 

DATE  AND  TIME  OF  BATHYTHERMOGRAPH  TRACE: 

DATE .  .10/23/63 

TIME .100832 

BATHYTHERMOGRAPH  SERIAL  NUMBER: 

BT  SERIAL  NO 663977 

NO  TEMPERATURE  INVERSIONS  EXIST  IN  THIS  TRACE. 

TRANSIENT  DATA  AS  COMPUTED  BY  PROGRAM: 

BEGINNING  OF  FIRST  TRANSIENT 21.0 

END  OF  TRANSIENT 55.0 

BEGINNING  OF  NEXT  TRANSIENT 60.0 

END  OF  TRANSIENT 87.2 

MULTIPLE  THERMOCLINE  DATA  AS  COMPUTED  BY  PROGRAM: 

CENTER  NUMBER  1  LOCATED  AT 148.0 

CENTER  OF  THERMOCLINE  149.6 

TOP  OF  THERMOCLINE 99.3 

BOTTOM  OF  THERMOCLINE  196.8 


END  OF  DATA  FOR  THIS  BATHYTHERMOGRAPH  SLIDE. 

NUMBER  OF  DATA  PAIRS  ANALYZED  -  18 


Fig.  8   Example  of  Output  for  Real  Data  Curve 
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START  OF  INVESTIGATIONS  IN  BATHYTHERMOGRAM  TRACE  NUMBER  27 

GEOGRAPHIC  POSITION  OF  BATHYTHERMOGRAPH  TRACE: 

LATITUDE   .....     00.00 

LONGITUDE 00.00 

DATE  AND  TIME  OF  BATHYTHERMOGRAPH  TRACE: 

DATE 00/00/00 

TIME 000000 

BATHYTHERMOGRAPH  SERIAL  NUMBER: 

BT  SERIAL  NO 000000 

THE  FOLLOWING  INDICATES  THE  COMPUTED  TOP  OF  THE  TEMPERATURE 
INVERSION,  BY  NON-GAUSSIAN  THEODS : 

TEMPERATURE  INVERSION  TOP  .....     66.0 

THE  FOLLOWING  INDICATES  THE  COMPUTED  BASE  OF  THE  TEMPERATURE 
INVERSION,  BY  NON-GAUSSIAN  METHODS: 

TEMPERATURE  INVERSION  BASE 160.1 

THE  FOLLOWING  FIGURES  DEFINE  THE  MINIMUM  TEMPERATURE  OF  THE 
TEMPERATURE  INVERSION: 

MINIMUM  TEMPERATURE  11.0 

DEPTH  OF  MINIMUM  TEMPERATURE 139.4 

THE  FOLLOWING  FIGURES  DENOTE  THE  TEMPERATURE  AND  SPACIAL 
MAGNITUDE  OF  THE  TEMPERATURE  INVERSION: 

THERMAL  MAGNITUDE   4.0 

SPACIAL  MAGNITUDE 94.1 

NO  THERMAL  TRANSIENTS  EXIST  IN  THIS  TRACE. 

NO  MULTIPLE  THERMOCLINES  EXIST  IN  THIS  TRACE. 

END  OF  DATA  FOR  THIS  BATHYTHERMOGRAPH  SLIDE. 
NUMBER  OF  DATA  PAIRS  ANALYZED  - 

Fig.  10  Example  of  Output  for  Inversion 
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Fig.11     Hand  drawn  Bathythemogram 
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START  OF  INVESTIGATIONS  IN  BATHYTHERMOGRAM  TRACE  NUMBER   17 
GEOGRAPHIC  POSITION  OF  BATHYTHERMOGRAPH  TRACE: 

LATITUDE   .....  

DATE  AND  TIME  OF  BATHYTHERMOGRAPH  TRACE: 

DATE   .....  

TIME 

BATHYTHERMOGRAPH  SERIAL  NUMBER: 

BT  SERIAL  NO. 

NO  TEMPERATURE  INVERSION  EXIST  IN  THIS  TRACE. 
NO  THERMAL  TRANSIENTS  EXIST  IN  THIS  TRACE. 
MULTIPLE  THERMOVLINE  DATA  AS  COMPUTED  BY  PROGRAM: 

CENTER  NUMBER  1  LOCATED  AT   ......... 

CENTER  NUMBER  2  LOCATED  AT   

CENTER  OF  GAUSSIAN  THERMOCLINE   

TOP  OF  GAUSSIAN  THERMOCLINE  ....  

BOTTOM  OF  GAUSSIAN  THERMOCLINE   

SEPERATION  #1  BETWEEN  THERMOCLINES   

CENTER  OF  GAUSSIAN  THERMOCLINE   

TOP  OF  GAUSSIAN  THERMOCLINE  

BOTTOM  OF  GAUSSIAN  THERMOCLINE  ..... 

END  OF  DATA  FOR  THIS  BATHYTHERMOGRAPH  SLIDE. 
NUMBER  OF  DATA  PAIRS  ANALYZED  - 


000.00  NORTH 
000.00  WEST 


00/00/00 
00000 


000000 


72.0 

144.0 

75.0 

53.5 

94.4 

100.0 

138.0 

104.9 

152.6 
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Fig.  12  Example  of  Output  for  Multiple  Thermoclines 
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COMPUTER    TRACE 
TEMPERATURE  AGAINST  DEPTH  (Z   POS.  UP) 

FIG.13 
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All  the  preceeding  traces  were  read  in  in  data  deck  format,  with 
varied  intervals  of  delta-T,  0.1  degrees  in  the  first  example,  1.0 
degrees  in  the  second  case,  and  0.5  degrees  in  the  latter  case.   There 
are  other  possible  combinations,  but  these  three  cases  outline  the  general 
format  and  the  answers  that  can  be  expected.   The  other  cases  are  just 
various  combinations  or  omissions  of  the  above  three  cases. 
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6.    Discussion  and  Conclusions  Concerning  Results. 

Although  no  error  analysis  was  undertaken  for  this  computer  pro- 
cess, it  is  obvious  that  there  are  some  errors  built  into  the  program 
that  could  perhaps  be  reduced  by  further  investigation  and  refinement. 
In  several  segments  we  use  a  Gaussian  based  method  and  thus  there  is 
some  error  introduced  in  the  computation  of  the  top  and  bottom  of  the 
thermocline  due  to  having  to  compute  a  sigma  in  the  program.   Also,  in 
this  program  the  multiplication  factor,  SK,  is  calculated  only  roughly 
based  on  the  thermal  interval.   This  factor  can  be  improved,  but  cannot 
be  exact  in  the  computer,  due  to  truncation  error;  besides,  much  storage 
space  is  needed  to  store  the  large  table  of  values  applicable  where  each 
tabular  figure  is  only  good  for  one  specific  Bet  of  data. 

There  are  certain  errors  inherent  in  the  computer,  but  since  we  do 
not  require  great  precision,  these  errors  can  be  ignored  due  to  the  fact 
that  the  bathythermograph  itself  is  not  a  precise  instrument.   These 
errors  include  truncation  and  round-off  errors  in  every  computation. 
Since  we  are  requiring  printout  to  only  one  decimal  place  it  is  doubt- 
ful that  these  errors  will  ever  be  noticable.  Also  there  are  errors 
built  into  the  functions  that  are  interral  to  the  computer,  such  as  the 
"square  root,"  "multiplication,"  and  "addition"  routines,  but  again  these 
will  not  be  significant  in  the  accuracy  that  we  require. 

There  are  certain  advantages  to  the  use  of  data  tapes  versus  magne- 
tic cards,  when  great  volumes  of  data  are  to  be  processed.   The  time 
required  to  punch  the  data  cards  and  to  check  the  results  is  a  function 
of  the  experience  of  the  operator.   The  computer  time  is  approximately 
25  seconds  for  card  input  computations,  plus  additional  time  if  the  draw- 
ing of  graphs  is  required  (approximately  three  minutes  per  graph).   This 
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time  can  be  cut  approximately  in  half  if  the  program  compiled  by  the 
computer  is  put  on  tape  and  data  input  is  in  that  form,  as  the  compiler 
takes  approximately  half  the  time  now  used  in  making  a  computer  run. 

There  are  other  equipment  errors,  that  as  yet  have  not  been  evalu- 
ated,that  are  concerned  with  the  input  data,  and  are  still  vital  to  the 
person  desiring  to  analyze  the  traces  by  use  of  this  program.   Such  items 
as  hysteresis,  where  a  double  trace  is  formed  on  the  bathythermograph 
slide,  due  to  slight  differences  in  the  temperature  when  raising  the 
Bathythermograph  as  compared  to  when  lowering  it  cause  inaccuracies, 
dull  styluses  can  give  traces  that  are  too  wide  to  be  useful.   Damaged 
equipment  that  gives  unacceptable  temperature  fluctuations,  such  as 
variations  from  bucket  thermometer  temperatures,  or  other  signs  that  the 
bathythermogram  is  inaccurate  although  the  trace  itself  appears  accept- 
able is  another  input  error.   Also  there  are  errors  inherent  in  the  digiti- 
zer that  can  vary  depending  upon  the  equipment  used  and  the  usual  problems 
inherent  to  electronic  curve  followers  or  other  devices. 

The  results  of  this  program,  in  the  fifty  test  cases  run,  appear  to 
give  slight  errors  as  compared  to  the  hand  analysis  of  the  same  data. 
Ten  of  the  test  cases  were  hand  constructed  curves,  designed  to  test 
the  program  for  specific  features  or  combinations  of  these  features,  and 
the  remaining  data  were  real  data  obtained  from  the  Fleet  Numerical  Weather 
Central  or  taken  from  other  written  works,  such  as  Tully  and  Dodimead 
(1957)  and  Pattullo  and  Cochrane  (1951),  and  a  few  traces  taken  in 
Monterey  Bay.   Both  tape  inputs  from  Fleet  Numerical  Weather  Central 
and  punched  tape  data  decks  were  used.   There  were  no  apparent  differ- 
ences between  the  outputs  from  the  digitized  data  and  the  data  on 
aperature  cards,  as  might  be  expected  assuming  that  care  was  taken  by 
the  person  picking  off  data  points  from  the  curves.   The  largest  error 
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in  the  tested  cases,  was  found  in  a  very  steep  thermocline,  real  data 
case,  where  the  depth  error  was  4.5  feet  for  the  top  of  the  thermocline 
and  3.0  feet  for  the  base  of  the  thermocline  in  a  thermocline  that  was  76 
feet  in  vertical  magnitude.   In  this  case  the  center  of  the  thermocline 
was  only  0.3  feet  from  the  hand  calculated  value.   The  thermal  interval 
was  1.0  degree  Farenheit  and  the  thermocline  in  question  was  analyzed 
by  the  Gaussian  technique. 

Thus  it  can  be  said  that  this  program  accomplishes,  with  acceptable 
accuracy,  considering  the  analytic  methods  used,  in  seconds  what  a  well 
trained  operator  would  take  tens  of  minutes  to  accomplish. 

This  program  can  be  expanded,  by  the  inclusion  of  additional  sub- 
routines, and  their  respective  "call"  cards  in  the  MAIN  routine  without 
greatly  increasing  the  time  needed  in  computing  the  desired  results. 
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7.    Recommendations . 

As  previously  stated,  there  are  other  features  that  can  be  added 
to  this  program  by  merely  including  additional  subroutines.   Additionally, 
this  same  programming  technique  can  be  used,  with  a  few  basic  changes,  to 
analyze  other  oceanographic  data  curves.   Before  accomplishing  any  con- 
versions of  this  program  to  handle  other  data,  the  user  must  first  look 
at  the  basic  form  of  the  traces  he  plans  to  analyze.  Traces,  other  than 
temperature,  may  not  be  describable  using  the  same  parameters  as  des- 
cribed by  Tully  (1964).   Such  items  as  the  magnitude  of  the  fluctuations 
of  other  quantities  must  be  examined  as  they  may  exceed  the  number  that 
this  program  can  handle.   For  example,  it  may  be  common,  in  some  other 
oceanographic  parameter,  to  have  more  than  three  stable  features  corres- 
ponding to  the  thermoclines  in  this  program,  and  thus  the  program  would 
have  to  be  modified  to  allow  for  additional  curve  computations.   Another 
correction  that  would  have  to  be  made  would  be  the  rewording  of  the  print- 
out statements. 

Further  investigations  should  be  undertaken,  to  include  the  improve- 
ment in  accuracy  of  sigma  and  SK  parameters  and  the  possibility  that  such 
curves  as  the  hyperbolic  tangent  or  Chi-squared  curves  might  give  a  better 
curve  fit. 

As  this  program  has  yet  to  be  run  on  the  CDC  6500  computer  of  the 
Fleet  Numerical  Weather  Central,  this  should  be  undertaken  as  soon  as 
possible.   For  this  type  of  program,  the  speed  of  computation  on  the  CDC 
6500  computer  is  approximately  three  times  that  of  the  IBM  360.   This 
CDC  6500  computer  also  contains  more  storage  space  than  the  IBM  360 
computer  and  thus  longer  tapes  or  card  inputs  may  be  used.   An  additional 
advantage  of  this  increased  storage  is  that  investigation  could  be  shifted 
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from  the  present  two  dimensional  set  up  (temperature  versus  depth)  to 

a  three  dimensional  investigation  (multiple  traces  analyzed  simultaneously) 

giving  a  field  representation  instead  of  just  specific  points. 
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Initialize  NBT 


* 

PRINT  OUT 

Print 

out  "Start  of  investigations  into 

trace 

characteristics  for  bathythermo- 

graph 

trace  number" 

[  CALL  SEE13    J 


READ:T(I),Z(I) 


YES 


Does  T(I)  equal  zero  and  simultaneously 
Z(I)  greater  than  zero? 


I 


Does  T(I)  equal  zero  and  simultaneously 
7,(1)    equal  zero? 


<•> 


CE 


I 


crement  counter 


I 


READ: 

LAT,XLAT,  LONG,  YLONG,  QUAD1,  0UAD2 , 
DATE1,  DATE2,  DATE3,  TIME,  SERIAL,  to 
determine  the  time,  serial  number,  and 
location  of  bathythermograph  drop 


Call  PRINT7 


] — •© 


NO 


(2> 


Test  for  repeat  values  of  temperature 
in  individual  trace. 


I 


Call   INVERS 


® 
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] — ■© 


(2> 


NO 


NO 


®" 


<8> 


f 


PRINT  OUT 

Print  out,  "No  thermal  transients  exist 
in  this  trace." 


I 


PRINT  OUT 

Print  out,  "No  multiple  thermoc lines 
exist  in  this  trace." 


I 


Test  thermal  interval  to  see  whether  or  not 
it  is  worthwhile  investigating  for  the 
presence  of  transients, 


I 


CALL  TRANS 


<5> 


Investigate  the  data  set  to  see  whether  there 
exist  more  than  one  major  thermocline.   If 
more  than  one  thermocline  exists,  each 
thermocline  is  investigated  as  to  the  location 
of  its  center,  top,  and  base.   If  only  one 
thermocline  exists  then  a  more  complete  analy* 
sis  is  made  of  this  single  thermocline. 


■•■ 


I 


Call   MANYT 


<■) 


PRINT  OUT 
Print  out,  "Completion  of  all  investiga- 
tions into  bathythermograph  trace  number." 


I 


55 


i 


Is  ICOUNT  greater  0 


\  END/ 


56 


L 


Subroutine    PRINT7 


J 


Inputs  from:  MAIN 
LAT,  XLAT,  LONG,  YLONG,  QUAD1,  0UAD2 , 
DATE1,  DATE2,  DATE3 ,  TIME,  SERIAL 


,YES 


I 


QUAD1-1 


3 


QUAD2-1 

ZZT" 


NO 


NO 


QUAD2-1 


j 


3= 


PRINT  OUT 

Print  out;  "NORTH*  latitude  and  "EAST" 
longitude  and  also  date,  time,  and  serial 
of  bathythermograph. 


PRINT  OUT 

Print  out;  "NORTH"  latitude  and  "WEST" 
longitude  and  also  date,  time,  and  serial 
of  bathythermograph. 


PRINT  OUT 

Print  out;  "SOUTH"  latitude  and  "EAST" 
longitude  and  also  date,  time,  and  serial 
of  bathythermograph. 


i 

® 
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PRINT  OUT 

Print  out;  "SOUTH"  latitude  and  "WEST" 
longitude  and  also  date,  time,  and  serial 
of  bathythermograph. 


RETURN 


5 
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NO 


Subroutine    INVERS 


I 


Inputs  from:  MAIN 
Z(I),  T(I),  ICOUNT 


I 


Define    SK   in    terms    of    temperature    interval 
Define    and    set   do-loop   counters. 


Call  0RDER2 


I 


=® 


Test  entire  data  set  to  Insure  inversion 
exists  by  checking  sign  changes  of: 
T(N/1)  -  (T(N)). 


I 


By  do-loop,  compute  first  difference,  A(N) , 
on  Temperature,  and  first  difference,  B(N), 
on  depth,  terminating  loop  if  sign  of 
A(N)  changes  ^^ 


By  do-loop,  compute  second  difference, 
G(L) ,  based  on  B(N). 


I 


Set  up  initial  testing  loop  and  set 
TEST1  equal  to  G(L) .   Run  loop  until 
sign  of  TEST1  changes  and  note  loop 
counter  at  this  poinj 


T 


Compute  top  of  inversion,  ZC1,  based  on 

TEST1  termination  counter,  and  also  comp- 

ute temperature  at  top  of  inversion, 

TEMP2. 

5 
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By  do-loop,  compute  first  difference, 
C(I31),  on  temperature,  and  first  differ- 
ence, D(I31),  on  depth,  terminating  loop 
if  sign  of  C(I31)  changes. 


T 


Call  ORDER  1 


I 


=© 


By  do-loop,  compute  second  difference, 
P(I32),  based  on  D(I31). 


I 


Compute  spacial  magnitude  of  inversion, 
XUPP,  based  on  loop  counter  NNCO. 


By  do-loop,  compute  first  difference, 
E(I33) ,  on  temperature,  and  first  differ- 
ence, F(I33),  on  depth,  terminating  loop 
at  counter,  based  on  length  of  data  set, 
ICOUNT. 


I 


Bv  do-loop,  compute  second  difference, 
R(I34)  on  F(I33). 


I 


Set  up  second  testing  loop  and  set 
TEST2  equal  to  R(I35).   Run  loop  until 
sign  of  TEST2  changes  and  note  loop 
counter  at  this  point. 


I 


Compute  base  of  inversion,  ZC2,  based 
on  counter  when  TEST2  sign  change 
encountered . 


60 


Compute  minimum  temperature,  TMIN,  and 
its  corresponding  depth,  ZMIN,  based  on 
counter  in  R(I34)  do-loop. 


I 


Call  PRINT3 


I 


Call  PRINT4 


I 


Call  PRINT5 


I 


Call  PRINT6 


I 


RETURN 
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(•) 


Subroutine  0RDER2 


I 


Inputs  from:  INVERS 
X(I),  III,  M 


I 


This  subroutine  takes  the  input  data  and 
rearranges  it,  if  necessary,  in  correct 
descending  order,  by  using  a  double  do- 
loop. 


I 


DO   I  -  M,  III-3 


I 


IP2  -  1/  1 


DO  I  -  IP1,  ni-3 
X(I).LE.X(J) 


SUBST   -  X(I) 

1 

X(I)    -     X(J) 

1 

X(J)      -   SUBST 

I 

CONTINUE 

I 


RETURN  H*- 


d2 


Subroutin 


e  0RDER1 


I 


I] 


Inputs  from:  INVERS 
X(I),  T(I),  III 


I 


This  subroutine,  similar  to  0RDER2,  takes  the 
input  depths  and  temperatures  and  sets  them  in 
correct  order  using  a  double  do-loop. 


I 


«■»    DO   I  -  1,  III-3 


I 


L 


IP1  -  1/1 


< 


I 


DO  J-IP1,  III-3 


:i 


X(I).LE.X(J) 


]_YES 

I  1 


SUBST  =  X(I) 


I 


1 


X(I)  -  X(J) 


I 


X(J)  »  SUBT 


SUBST  «  T(I) 


1 


b3 


T(I)  -  T(J) 


T(J)  -  SUBST 


CONTINUE 


RETURN 
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Subroutine    PRINT3 


I 


Inputs  from:   INVERS 
ZC1 


I 


PRINT   OUT 
Print    out,    "top   of    temperature    inversion 
based    on   on-gaussian  methods" 


I 


RETURN 
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Subroutine      PRINT4 


I 


Inputs  from:   INVERS 
ZC2 


PRINT  OUT 
Print  out,  "Base  of  temperature  inversion 

based  on  non-gaussian  methods." 


I 


RETURN 


66 


Subroutine    PRINT5 


I 


J 


|  1 

Inputs    from:       INVERS 

'  TMIN ,    ZMIN 


I 


PRINT  OUT 

Print  out,  "Minimum  temperature  in  tempera- 

ture inversion"  and  "Minimum  temperature 

found  at  depth." 

— 

I 


RETURN 


1 
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0 


Subroutine    PRINT6 


I 


Inputs  from:  INVERS 
ZUPP,  TUPP 


I 


Print  out,  "Spacial  magnitude  of  temperature 
inversion"  and  "Thermal  magnitude  of  tempera- 
ture inversion." 


I 


RETURN 
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Subroutine   TRANS 


I 


Inputs  from:   MAIN 
Z(I),  T(I),  ICOUNT 


I 


Set  do-loop  counters,  based  on  size  of 
data  set,  ICOUNT. 


i 


Bydo-loop   compute    first   difference, 
Z6(N) ,    on  depth. 


I 


By  do-loop,  compute  second  difference, 
Z7(I^1) ,  on  Zb(N). 


I 


.YES  _ 


By  loop,  check  sign  of  second  difference 
to  see  if  it  is  greater  than  zero. 


I 


PRINT  OUT 
Print  out,  "Beginning  of  first  transient" 


YES 


I 


NN.LE.M 


I 


CONTINUE 


NN.LE.M 


YES 
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YES 


YES 


-•- 


Increment   NN 


I 


Z7(NN).gt.O? 


Increment  NN 


I 


NN.LE.M? 


I 


Increment  NN 


i 


] 


Z7(NN).LT.O? 


I 


INCREMENT  NN 


NN.LE.M? 


I 


PRINT  OUT 
Print  out,  "End  of  transient." 


I 


YES 
» 


YES 


PRINT  OUT 
Print  out,  "Beginning  of  next  transient" 


I 


NN.LE.M? 


70 


® 


RETURN 


71 


£ 


Subroutine  MANYT 


I 


Inputs  from:   MAIN 
Z(I),  T(I),  ICOUNT 


I 


I 


Initialize  counters 


Set  do-loop  counters  based  on  size  of 
data  set  as  inputed. 

1 

1 

By  do-loop,  compute  first  difference, 
Z6(N2),  on  depth. 

I 

By  do-loop,  compute  second  difference, 
Z8(N1),  on  Z6(N2). 

By  do-loop,  check  all  values  of  Z8(N1) 
to  see  if  greater  than  zero. 


Increment  cou 


nter     j 


PRINT  OUT 
Print  out  "Center  number  , 


located  at' 


LZZ 


N.LT.M? 
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® 


Set  value  to 
BIG1 


By  do- loop,  compute  the  largest  value  of 
Z6(n)  in  therraocllne.   Set  BIGDP1  equal  to 
0.5  (Z6(max)  +  Z6(max  +1)). 


I 


Call  GAUSS3 


<2> 


PRINT  OUT 
Print  out,  "separation  #1  between  thermo- 
clines." 


I 


Reinitialize  counters 


I 


Set  value  to: 
BIG2 


I 


By  do-loop,  compute  the  largest  value  of 
Zb(N)  in  the  next  thermocline.   Set 

BIGDP1  equal  to  0.5(Z6(max)  +  Z6(max  +  1)) 


1 

Call  GAUSS3 

■» 

^           1 

1 

1 

PRINT  OUT 
Print  out,  "Seperation  #2  between  thermo- 
clines." 

5 
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YES 


Rein 


itialize  counters 


I 


D 


CALL  GAUSS3 


I 


=® 


If  number  of  the rmoc lines  is  greater  than 

one . 


I 


If  there  is  only  one  thermocline  in  the  data 
set,  then  call  NONGA. 


I 


P® 


PRINT  OUT 

Print  out,  "End  of  data  for  this  bathy thermo- 
gram trace"  and  "Number  of  data  pairs 
analyzed .  " 


I 


RETURN 


£ 
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Subroutine  GAUSS3 


Inputs  from:   MANYT 
Z(I),   T(I),   ICOUNT,  NQN1 


Compute  value  of  SK,  based  on  the  tempera- 
ture interval. 


I 


Call  0RDER3 


I 


t=® 


Compute  Z1(N)  -  (Z(N)  +  Z(n+l)/2.0)  and 
Z2(N)  -  Z1(N)  x  Z1(N). 

V 

Sum:   Z1(N)  and  Z2(N)  within  do-loop. 

1 

f 

Compute  gaussian  center  of  thermocline, 
ZC  -  (1.0/(total  point  -  l))xsum  Z1(N). 

< 

! 

PRINT  OUT 
Print  out  "Center  of  Gaussian  thermocline 
is  located  at  depth." 

' 

1 

Compute:   ZAVG  ■  (1.0/(total  -1)  x  sum 
of  Z2(N);   SIGMA  =■  ZAVG  -  ZCxZC; 

I 


Compute  top  of  Gaussian  thermocline 
ZT  -  ZC  -  SK  x  SIGMA 


S 
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® 


PRINT  OUT 


Print  out  "Top  of  gaussian  thermocline  is 
located  at  depth" 


Compute  base  of  gaussian  thermocline, 
ZBT  »  ZC  +  SK  X  SIGMA 


1 


PRINT  OUT 
Print  out  "Base  of  gaussian  thermocline  is 
located  at  depth" 


I 


Compute  magnitude   of    thermocline,    from 
WW   *   ZBT   -    ZT. 


I 


PRINT  OUT 
Print  out  "Magnitude  of  thermocline." 


RETURN 
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s 


Subroutine  0RDER3 


I 


Input 8  from:   NONGA 
Z(I),  KT 


I 


This  subroutine  arranges  the  input  data, 
if  necessary  into  correct  descending  order 
from  the  surface  down. 


I 


DO   I  -  1,  KT-1 


I 


IP1  -  I  +  1 


3=Z 


DO   J  -  IP1.  KT-1 
~j 


X(I).LE.X(J) 


I 


SUBST  -  X(I) 


I 


3 


X(I)  -  X(J) 


I 


X(J)  -  SUBST 


I 


CONTINUE 
RETURN 


YES 
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Subroutine  NONGA 


Inputs  from:   MANYT 
Z(I),   T(I),    ICOUNT 


I 


Compute  SK ,  based  on  the  temperature 
inte  rval . 


I 


Compute  first  difference,  by  do-loop, 
Z1(N) ,  on  depth. 


I 


Compute  by  do-loop,  second  difference, 
Z2(n),  on  Z1(N). 


] 


i 


Set  up  testing  loop,  TEST(KK)  equal  to 
Z2(KK)  and  execute  loop  until  there  is 
a  change  in  sign  of  TEST  (KK) . 


1 


Compute  the  center  of  the  thermocline, 
ZC,  based  on  the  counter  number  when 
TEST  loop  terminated. 


I 


Set  up  localization  loop  to  find  extent  of 
the  upper  portion  of  the  thermocline 


Call  GAUSS1 


I 


J= 


Set  up  localization  loop  to  find  extent  of 
the  lower  portion  of  the  thermocline 
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Call      GAUSS2 


r 


*=@ 


RETURN 
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I        Subroutine   GAUSS  1 


I 


Inputs  from:   NONGA 
Z(D,   NT,    ICOUNT 


I 


] 


Compute  value  of  SK,  based  on  input 
value  of  ICOUNT. 


L 


J 


Compute;  Z1(N)  =  (Z(N)  +  Z(N+l))/2.0 


I 


Compute;  Z2(N)  »  Zl(n)  x  Z1(N) 


I 


3 


Increment  N 


I 


By  do-loop,  compute  SUM1  and  SUM2,  based 
on  sums  of  Z1(N)  and  Z2(N). 
I , 1 


I 


Compute;    ZC   »    (1.0/total    -1)    x    (SUM1),    the 
center   of   the   gaussian   thermocline. 


Compute;    ZAVG   «    (1.0/total    -    l)x(SUM2); 
SIGSQR   -   ZAVG    -    (ZCxZC) ;    ZT  =   ZC  +  SIGSQR 


J 


Call    PRINT1 


i 


fc^@ 


RETURN 
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Subroutine  GAUSS 2 


Inputs  from:   NONGA 
Z(I),   NT,   ICOUNT 


Compute  value  of  SK,  based  on  input  value 
of  ICOUNT. 


Compute   Z1(N)    -    (Z(N)    +  Z (N+l) ) /2.0 


Compute    Z2(N)    -    (Z1(N)    x   Z1(N)) 


I 


Increment   N 


I 


i 


Compute;  ZAVG  =»  (1.0/total  -1)  x  SUM2 ; 
SIGSQR  -  ZAVG  +   (ZC  x  AC). 


I 


Compute  base  of  thermocline,  ZBT  -  ZC  + 
SK  x  SIGSQR.  


By   do-loop,    compute   SUM1   and    SUM2   based 
on    sumes    of   Z1(N)    and    Z2(N). 

■  -    _ 

I 

' 

Compute;      ZC»    (1.0/total    -    1)    x   SUM1 

Call   PRINT2 


=Tae 
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Subroutine  PRINT1 

Input 8  from;  GAUSS  1  1 
ZT,  ZC,  ICOUNT 


PRINT  OUT 
"Number  of  data  pairs  utilized" 


I 


PRINT  OUT 

Print  out  "Center  of  thermocline  is 
located  at  depth. 


I 


PRINT  OUT 

Print  out  "Top  of  thermocline  is 
located  at  depth." 


3g 

RETURN  | 
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Subroutine   PRINT2 


1_J 


i 


Inputs    from:      GAUSS2 


ZBT 


L_ 

PRINT  OUT 


Print  out  "Base  of  thermocline  located  at 
depth." 


1 


!  RETURN 


: 
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APPENDIX  II 
OBJECTIVE  THERMOCLINE  ANALYSIS  PROGRAM 


84 
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APPENDIX  III 
WRITE-UP  OF  "DRAW"  SUBROUTINE 
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